Disorder induced quantized conductance with fractional value and universal 
conductance fluctuation in three-dimensional topological insulators 
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We report a theoretical investigation on the conductance and its fluctuation of three-dimensional 
topological insulators (3D TI) in Bi2Sez and Sb2Te3 in the presence of disorders. Extensive numer- 
ical simulations are carried out. We find that in the diffusive regime the conductance is quantized 
with fractional value. Importantly, the conductance fluctuation is also quantized with a universal 
value. For 3D TI connected by two terminals, three independent conductances Gzz, Gxx and Gzx 
are identified where z is the normal direction of quintuple layer of 3D TI (see inset of Fig.l). The 
quantized conductance are found to be (Gzz) = 1, {Gxx) = 4/3 and (Gzx) = 6/5 with corresponding 
quantized conductance fluctuation 0.54, 0.47, and 0.50. The quantization of average conductance 
and its fluctuation can be understood by theory of mode mixing. The experimental realization that 
can observe the quantization of average conductance is discussed. 

PACS numbers: 72.10.Fk, 73.20.Fz, 73.25. -(-i, 



CD 

S 

a 

o 
o 



> 
cn 

00 
00 

m 
d 



X 



Recently, the topological insulator (TI), a new state 
of matter, has attracted a lot of theoretical and experi- 
mental attention.—"— The TI has an insulating energy gap 
in the bulk states which behaves like the general insula- 
tor, but it has exotic gapless metallic states on its edges 
or surfaces. The TI is first predicted in two-dimensional 
(2D) systems, e.g., the graphene and HgTe/CdTe quan- 
tum well. It has been generalized* in 3D and confirmed 
experimentally.— The 2D TI has the gapless helical edge 
states and exhibits the quantum spin Hall effect while in 
3D TI the conducting state is helical surface state. This 
hehcal edge or surface states are topologically protected 
and are robust against all time-reversal-invariant impu- 
rities. Many interesting physical phenomena have been 
predicted including Majorana fermion^, topological mag- 
netoelectric effect^, magneto-optical Kerr and Faraday 
effects.— There are also many studies on disordered TI. It 
was found by Li et al that in the presence of disorder,^ 
a new phase called topological Anderson insulator can 
be induced. The complete physical picture and mecha- 
nism of topological Anderson insulator (TAI) was given 
by Groth et al^'^ using an effective medium theory. Re- 
cently, the TAI was also predicted in 3D TlJ^ Both ana- 
lytic and numerical results show that disorder can induce 
strong topological Anderson insulator in 3D. 

Disorder can affect mesoscopic systems in an important 
way. For instance, it is well known that, in the diffusive 
regime, the conductance fluctuation of the mesoscopic 
system assumes a universal value that is independent of 
system parameters and depends only on dimensionality 
and symmetry of the system For a 2D TI (HgTe/CdTe 
quantum well), it was found that the spin-Hall conduc- 
tance fluctuation is universal. ^'^ It will be of great interest 
to explore the universal behavior of conductance and its 
fluctuation in 3D TI. 

In this paper, we have studied the effect of disorder 
on topological surface states in 3D TI for Bi^Se^ system 
using extensive numerical simulation. Due to anisotropy 
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FIG. 1: (Color online) Conductance and its fluctuation vs 
disorder strength for different transport directions, Gzz for 
panel (a) and (6), Gxx, Gyy and Gxy for panel (c) and (d), 
Gxz for panel (e) and (/). Each data point on the flgure is 
averaged over 5000 configurations. 



of the Bi2Se3 system, there are three independent two- 
terminal conductance G^^, Gxx and Gzx where z is the 
normal direction of quintuple layer of 3D TI. We found 
that in the diffusive regime, the average conductance is 
quantized with fractional value 1,4/3,6/5, respectively, 
for Gzz, Gxx and Gzx (in units of e^/h). The correspond- 
ing conductance fluctuation is also quantized. By vary- 
ing the disorder strength W and the Fermi energy Ep 
inside the bulk gap, the phase diagrams of conductance 
and its fluctuation in the plane of {W,Ef) were con- 
structed, showing the same quantized conductance and 
fluctuation for a wide range of Ep and W . By introduc- 
ing the number of effective transmission and reflection 
channels, the quantized conductance as well as the quan- 
tized fluctuation can be understood by the mode mixing 
theoryJ^ The quantized values of conductance and fluc- 
tuation from the proposed effective mode mixing theory 
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are in excellent agreement with our numerical results. It 
is known that as disorder increases the electron transport 
goes through the ballistic and the diffusive regime and fi- 
nally enters the localized regime. Our numerical result 
shows that the evolution of topological surface states in 
the presence of disorders undergoes three stages: (1). the 
topological surface states are protected up to a critical 
disorder strength beyond which the topological surface 
states are destroyed while the bulk gap is present. At 
this stage, the conducting channels are non-topological 
surface states. This is true when the system begin enter- 
ing the diffusive regime. (2). While still in the diffusive 
regime but at large disorders, the bulk gap is destroyed 
and the conducting channels are bulk states. (3). For 
large enough disorders, the system becomes localized and 
there is no conducting channel available. Furthermore, 
four-terminal conductance and its fluctuation were cal- 
culated. Similar quantization behaviors for conductance 
and fluctuation were found. 

In our numerical simulation, we discretize spatial coor- 
dinates of the continuous effective k space Hamiltonian 
Hti in RefUl on the square lattice: 

Hti = J2 + J2 + H.c.il) 

where i = {ix,iy,iz) is the site index and 5x,5y,6z 
are unit vectors along x, y and z directions where z 
is the normal direction of quintuple layer of 3D TI. 

= (oi, 6i, Ci, di)"^, and ai,6i,Ci,di represents the four 
annihilation operators of electron on the site i with 
the state indices \P\+ ,\),\P2- ,\),\P\+ ,i),\P2- ,\). In 
Eq. ([T]) , Hii and H.^ are 4x4 Hamiltonian that 

are given by Ha = (C + 2Di/a? + AD2/ 



(M - 2Bi/a2 _ 4S2/a2)ro, H. g, 
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Sy 
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B2/a^To-iA2/i2a)Ti,H.j^ = -Dja^h^i + Bi/a^To- 

i^i/(2a)r2, where Foa,2 = (/2x2(8)s^, Sy(g)/|x2: Sz'8'/fx2) 
and is anti-diagonal identity matrix. Static Anderson 
type disorder is added to the on-site energy with a uni- 
form distribution in the interval [—W/2,W/2] where W 
characterizes the strength of the disorder. Here a is the 
lattice constant, and Ai, A2, Bi, B2, CJJi , D2 and M are 
system's parameters taken from Ref ij. 

By using Green's functions, the charge conductance 
from the terminal-/? to the terminal-a can be calcu- 
lated by using Landauer-Biittiker formula Gap{E) = 
{e^ / h)Tap , and = TrpaG'T^G"] is the transmission 
coefficient. The linewidth function TaiE) — i[E]^ — S°] 
and the Green's functions G'^^°'{E) can be calculated 
from C = [6"^]+ = [EI -Hc-J^^ S^-^, where He is 
Hamiltonian matrix of the central scattering region and / 
is the unit matrix with the same dimension as that of He, 
are self energy of external leads and can be calculated 
numerically''^^. The conductance fluctuation is defined as 
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rms(G) = Y (G^) — (G) , where (• • •) denotes averaging 
over an ensemble of samples with different disorder con- 



FIG. 2: (Color online)Averaged DOS in the middle slices of 
the central simulation box with injected wave from left lead 
with different disorder strength. 



figurations of the same strength W. In the following the 
average conductance and its fluctuation are measured in 
unit oi e^/h. 

In the numerical calculations, we choose the realistic 
material parameters for Bi2Se3^^ and perform calcula- 
tions on a L X L X L cubic sample with L = 20 and 
the lattice constant a = 5 A. Since each site has four 
orbitals, the dimension of the Hamiltonian or Green's 
function becomes 32000 for this system making the simu- 
lation very computational demanding. We first study the 
two terminal device by considering different transport di- 
rections in the presence of Anderson type disorder. As 
expected, our numerical simulation shows that transport 
properties along x direction is the same as that of y di- 
rection. So for the two-terminal structure there are four 
conductances: Gzz, Gxx, G^y, and Gzx that are plotted 
in Fig.l along with their fluctuations against the disorder 
strength. Here the Fermi energy is inside the bulk gap so 
that only topological surface states are conducting chan- 
nels. In general, our results show that the topological 
surface states are gradually destroyed in the presence of 
small disorders. As the strength of disorder increases, the 
system enters the diffusive regime and the conductance 
becomes quantized for a wide range of disorder strength. 
Importantly, this quantization of conductance is accom- 
panied by the quantized conductance fluctuation. This 
conductance fluctuation is universal since it is indepen- 
dent of parameters such as Fermi energy and disorder 
strength as shown in Fig. 3. From Fig.l, we see that 
the quantized conductance takes fractional value with 
{Gzz) = 1, (Gxx) = (Gxy) = 4/3, and (C) = 6/5. 
Two points worth mentioning from Fig. 1(c) and (d): (1). 
Two set of curves for (G) and rms(G) versus disorder 
strength along x direction and y direction are exactly the 
same which is expected. (2). For Gxx and Gxy, however, 
only their average conductance and its fluctuation in the 
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FIG. 3: (Color online) Phase diagram of conductance (panel 
(a)) and its fluctuation (panel (b)) in the plane of {W,Ef)- 
1000 samples are collected in each point on the phase diagram 
calculation. 



diffusive regime are the same. Its origin can be under- 
stood using the theory of mode mixing discussed below. 
To make sure that the quantization plateau lies in the dif- 
fusive regime, we have calculated the localization length 
in the quantization plateau region for W = [5, 10]. Our 
results give ^/L = 6~8in this range of W indicating 
that the system is indeed in the diffusive regime. 

The quantization of conductance in the diffusive 
regime can be described using the following phenomeno- 
logical theory of mode mixing in the central region, For 
instance, for the transport along z direction, two incom- 
ing topological surface states are completely mixed in 
the scattering region in the diffusive regime with equal 
probability of going forward pf and backward pb, i.e., 
Pf = Pb = 1/2. Hence the conductance is {Gzz) — '^Pf ~ 
1. Alternatively, we can use the number of effective 
transmission channels vt and reflection channels Vr with 
t^t/i^r — Pf/Pb- For Gzz, we have i^r ^ i^f = I so that 
Pf — vt/{i^t + i^r) = 1/2. For Gxx, only partial mixing is 
achieved with 1/^ = 1 and vt = 2, or pf ~ IjT). As a result, 
{Gxx) = 2pf = 4/3. Since x and y directions are equiva- 
lent, an incoming electron from x direction traversing in 
the scattering region find the same available vt along y di- 
rection as that of X direction. Hence (Gxx) = (Gxy)- To 
support this argument, we have calculated two-terminal 
conductance of a four-terminal device with four terminal 
on x-y plane (see the inset of Fig. 4(b)). For an electron 
from L-lead, we should have (Grl) — {Gul) — {Gdl) if 
x and y directions are equivalent. This is indeed what we 
found numerically. In addition, the conductance fluctu- 
ations are also the same. For Gzx, Vr — ^ and vt = 3/2 
due to the partial mixing, or (Gzx) = 2pf = 6/5. In 
terms of the number of effective channels, the quantized 
conductance can be given by the following expression. 
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FIG. 4; (Color online) Conductance and its fluctuation vs 
disorder strength for difi^erent transport direction in four ter- 
minal cross-setup. 



This ansatz gives the average conductance that agrees 
with our numerical result for the quantized conduc- 
tance of two-terminal structures. The universal conduc- 
tance fluctuations(UCF) can also be expressed in terms 
of the number of effective transmission and reflection 
channels 



rms(G) = 4< 



,2,, 2 



(3) 



(G) 



, VtVr 
Vt + Vr 



(2) 



With {ur,yf) = (1,1), (1,2), and (1,1.5) for Gzz, 
Gxx, and Gzx, respectively, we find from Eq.(j3l) that 
rms(G^^) = 2/VT4 = 0.535, rms(G^^) = 8/3/v^ = 
0.457, and rms(G^^) = 12/5/^23 = 0.500 that are very 
close to our numerical results (see Fig.l). 

We now examine the evolution of wave function as 
disorder strength is varied in the diffusive regime when 
the propagation is along z direction. In our numerical 
simulation we set Ep — 0.1348ey in Fig. 2(a), (b) and 
Ep — O.OlTSey in Fig. 2(c), (d). In the clean sample there 
are two pairs of topological surface states localized on the 
surface and separated in space so that the backscattering 
is forbidden. In Fig. 2(a), we have shown the DOS of one 
of the topological surface states propagating along posi- 
tive z direction. As the disorder is increased, these local- 
ized surface states become extended on the surface and 
the topological surface states are destroyed. However, 
current carrying states are still surface states and there is 
no DOS in the bulk (see Fig. 2(b)). This behavior of bulk 
insulator persists when the average conductance becomes 
quantized in the diffusive regime, e.g., when W = 6eV 
the system is a bulk insulator but the surface states is 
not topological protected. As disorder strength increases 
further the surface states gradually expand toward the 
center of the system and eventually become bulk states. 
At = 8eV, the bulk gap is closed due to the disor- 
ders. In order to investigate the nature of mode mixing. 
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we calculate the DOS for electron coming from the left 
lead along z direction at W — 8eV. Fig. 2(c) shows the 
average DOS in x — y plane where 10000 configurations 
are collected. From Fig. [51 we see that two modes are 
fully mixed in x — y plane due to the disorder scatter- 
ing. In the X — z plane (Fig. 2(d)), the DOS is mainly 
distributed around the left region due to the suppression 
of transmission ((G) ~ 1). 

To demonstrate the universal feature of the quantized 
conductance and its fluctuation, we have calculated the 
phase diagram where G^z and rms(Gz2) are plotted in 
{W,Ep) plane (see Fig. 3). In the calculation, the periodic 
condition in x direction is employed. The light green re- 
gion shows the quantized conductance. For Fermi energy 
inside the energy gap in the clean sample, we see that 
the system goes quickly from topological insulator to the 
region of quantized conductance. As disorder strength in- 
creases, it slowly enters the localized regime. This phase 
diagram shows that the conductance and its fluctuation 
are independent of Fermi energy and disorder strength in 
a wide range. 

Now we study the transmission coefficient T^^ of a 
four-terminal device shown schematically in the inset of 
Fig. 4(b) where the leads L and R are along z direction 
and leads U and D are along x direction. In the calcula- 
tion we have fixed Ep — COlTSeF and over 5000 samples 
are collected for each average. Fig. 4 shows that in the 
diffusive regime, the conductance and its fluctuation are 
again quantized. When there are four terminals, how- 
ever, the quantized conductance is no longer a fractional 
value. It is much smaller than that of the two-terminal 
case. This is not unexpected since there are two more 
terminals that electron can exit. To our surprise, the 
conductance fluctuation is also much smaller than that 
of the two-terminal case, nearly halved. 

To provide further evidence of universal conductance 
fluctuation, we have performed calculation for another 
3D TI in Sb2Te^ on a 20 x 20 x 20 cubic sample. Our 
numerical results give the same values of quantized con- 
ductance and fluctuation. So far, we have studied the 
influence of bulk disorder on the transport of 3D TI. We 



have also calculated the averaged conductance and its 
fluctuation in the presence of "surface" disorder. Specifi- 
cally, we have calculated Gzz when disorders are present 
only on the first layer or first three layers of surface of 3D 
TI in Bi2Se3. The quantized conductance and its fluctu- 
ation again show the same quantized plateau as the case 
of the bulk disorder with the same quantized conductance 
and quantized fluctuation. Note that the conductance 
measurement with surface disorder has been carried out 
experimentally on Bi2Se^!^— Since the quantization of 
conductance and fluctuation exist in a large window of 
disorder strength (corresponding to doping concentration 
in experiment), we believe that our results can be checked 
experimentally. 

To summarize, we have carried out extensive numer- 
ical simulation to calculate the transport properties of 
disordered 3D TI. Our results show that in the diffu- 
sive regime, the two-probe conductance is quantized with 
fractional value 1,4/3,6/5. The corresponding conduc- 
tance fluctuation is also quantized. An effective mode 
mixing theory is proposed that gives quantized conduc- 
tance and fluctuation in excellent agreement with our 
numerical results. The numerical results from different 
parameters including disorder strength and Fermi energy 
(from phase diagram of conductance and its fluctuation) , 
types of disorders (bulk and surface disorder), and types 
of 3D TI {Bi2Se3 and Sb2Te3) suggest that the quan- 
tized conductance and quantized fluctuation is a univer- 
sal property of 3D TI. We have also studied effect of dis- 
order on topological surface states. We found that as dis- 
order strength increases the conducting channels changes 
from topological surface states to non-topological surface 
states and finally to bulk states in the diffusive regime. 
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